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ABSTRACT 

We have analysed oscillations of the red giant star HD 186355 observed by the 
NASA Kepler satellite. The data consist of the first five quarters of science op- 
erations of Kepler, which cover about 13 months. The high-precision time-series 
data allow us to accurately extract the oscillation frequencies from the power 
spectrum. We find the frequency of the maximum oscillation power, Umax, and 
the mean large frequency separation, Av, are around 106 and 9.4 //Hz respec- 
tively. A regular pattern of radial and non-radial oscillation modes is identified by 
stacking the power spectra in an echelle diagram. We use the scaling relations of 
Ai/ and u^ax to estimate the preliminary asteroseismic mass, which is confirmed 
with the modelling result {M — 1.45 ± 0.05 Mq) using the Yale Rotating stellar 
Evolution Code (YREC7). In addition, we constrain the effective temperature, 
luminosity and radius from comparisons between observational constraints and 
models. A number of mixed I — 1 modes are also detected and taken into account 
in our model comparisons. We find a mean observational period spacing for these 
mixed modes of about 58 s, suggesting that this red giant branch star is in the 
shell hydrogen-burning phase. 

Subject headings: stars: individual (HD 186355) - stars: oscillations - stars: modelling 
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Introduction 



Stu dying solar-like oscillations p rovides a powerful method to probe the interiors 



of stars (jChristensen-Dalsgaard 



20041 ). Solar-like oscillations are expected in low- mass 



main-sequenc e stars cooler than the red edge of the classical instability strip in th e 



HR diagram 



Houdek et al. 



Christensen-Dalsgaard 



1982 



Christensen-Dalsgaard fc Frandsen 



1983 



19991). as well as in more evolved red giants which represent the future of 



our own Sun ( iDziembowski et al. 



2001 



Dupret et al. 



20091 ). It is thought that turbulent 



convective motions near the surface excite the oscillations stochastically. 



Asteroseismology of red giants has developed rapidly. It began with several detections 
of solar-like osc illations in G and K-type giants based on gro und-based observ ations in bot 



radial velocity (IFrandsen et al. 



2002 



De Ridder et al. 



20061 ) and photometry ( IStello et al. 



20071 ) and on space-based photornetry detections observed by the Hubble 



{HST- 



Edmonds fc Gilliland 



Infrared Explorer ( WIRE] 



1996 



Gilliland 



Buzasi et al. 



Mass Ejection I mager (SMEI; 



2000 



of Stars {MOST; iBarban et al 



arrant et al 



2007 



2008 



Stello fc GiUiland 



Retter et al. 



2003 



Space Telescope 



2009), Wide Field 



Stello et al 



20081 ) ■ Solar 



20071). and Mic rovariability and Oscillations 



Kallinger et al 



2008a 



b|). The oscillation periods 



in red giants range from hours to days. Ground-based observations usually suffer from 
interruptions and aliasing which complicate the measurement of oscillations. On the other 
hand, observations from space can provide high signal-to-noise ratio (SNR) and continuous 
data sets from which we may extract the oscillation parameters accurately. The 150-day 
observations by the Convection Rotation and planetary Transits sat ellite (CoRoT) clearly 



detected radial and non-radial oscillati ons in the range 10-100 /iHz (jPe Ridder et al. 



Hekker et al. 



2009 



Carrier et al. 



20091 : 



20101 ). which greatly increased the number of detected 



pulsating G and K giants and led to a huge breakthrough in the study of r ed giants. Thes e 



observations were followed by even more impressive results by Kepler (e.g. 



Bedding et al. 
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2010a: 


Huber et al. 


2010; 


Kallineer et al. 


2010; 


Hekker et al. 


2011a 


b; 


Di Mauro et al. 


2011 


Chaplin et al. 


2011 


)■ 



This paper presents observations and models of HD 186355 (HIP 96878, KIC 11618103), 
which is one of the brightest red giants in the Kepler field (V = 7.95). 



2. Observations 



The Kepler Mission (iBorucki et al. 



2008 



20101 ) was successfully launched on March 



7, 2009. Its primary scientific goal of is to search for Earth-sized planets in or near the 
habitable zone and to determine how many stars have this kind of planets in our Milky 
Way. Kepler is equipped with a 0.95-meter diameter telescope with an array of CCDs 
which continuously points to a large area of the sky in the constellations Cygnus and Lyra 
to detect the transits of the planets. Over the whole course of the mission (at least 3.5 
years), the spacecraft will measure the variations in the brightness of more than 100,000 
stars, which will be outstanding data for the study of asteroseismology. For many of 
these stars we can detect solar-like oscillations, which will allow us to investigate them in 
detail and obtain their fundamental properties, by us ing the techniques of asteroseismology 



(IChristensen-Dalsgaard et al. 



2007; 



Aerts et al 



20101 ) 



We used the first five quarters of data of HD 186355, which covers a total of about 
13 months. The raw long-cadence data (29.4 minutes sampling; Jenkins et al. 2010) 
were corrected by performing a point-to-point sigma clipping to remove the outliers. 
Additionally, a thermal drift was corrected by fitting a second-order p olynomial to the 



affected parts of the time-series. From the parallax of 5.44 ± 0.63 mas ( 



and using a bolometric correction for G5 giants of -0.34 from 



van Leeuwen 



20071 ) 



Kalerl (119891 ). we derived the 



luminosity of the star to be 24.0 ± 5.6 Lq. We take the effective temperature (Tgfj = 4867 
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± 150 K) from Kepler Input Catalogue (KIC; Brown et al. 2011). 



3. Global oscillation analysis 



Solar-like oscillations are usually high-order and low-degree p-mo des. Their frequencies 



1980 



Gough 



are r egularly spaced, approximately following the asymptotic relation (iTassoul 
19m : 

z/„; ^ Az/(n + i/ + e) - /(/ + l)Do, (1) 

where n is the radial order and / is the angular degree. The quantity Au (large frequency 
separation) is approximately the inverse of the sound travel time across the star, while e is 
sensitive to the surface layers and, for relatively unevolved stars, Dq is sensitive to the sound 
speed gradient near the core. As the star evolves, the stellar envelope starts to expand 
and the p-mode frequencies gradually decrease while oscillations in the core driven by 
buoyancy (g- modes) shift to higher frequencies. This eventually leads to so-called "mixed 
modes". These are non-radial oscillation modes that have a mixed character, behaving 
like g-modes in the core and p-modes i n the envelope, and shifting in freq uency as they 



undergo the so-called avoided crossings (jOsaki 



1975 



Aizenman et al. 



19771 ). For red giants. 



the asym ptotic / = 1 modes in particular depar t from the relation due to many avoided 
crossmgs fiHuber et allLoiOi : iMosser et al.l 



201l|). 



Our frequency analysis covers three basic steps that are performed on the power 
spectrum of the Kepler light curve: fitting and correcting for the background; estimating 
the frequency of maximum power (z/max) and the large separation (Az/); and extracting 
individual frequencies (uni)- In the following subsections, we describe the three analysis 
steps in detail. 
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3.1. Modelling the Background and Determining z/max 

The power spectrum shows a frequency- dependent background signal due to stellar 
activity, g ranulation and faculae which can be modelled by a sum of several Lorentzian-like 



functions (IHarvey 



19851 ). In this paper, stellar activity, granulation and faculae w ere 



KarofB (120081 ). which 



represented by modified Lorentzian-like functions, first introduced by 
give a better fit to the background than a Harvey model with a constant slope of —2. This 
background model has a shallower slope at low frequencies and a steeper slope at higher 
frequencies, corresponding to stellar activity and granulation, respectively. The power 
excess hump from stellar oscillations is approximately Gaussian, so the complete spectrum 
was modeled by: 

Pi^) =Pn + t ^tli, ,4 + ^.-P , (2) 

1 + [2-nvTiy + {2'nvTiY \ 2a^ J 

where P„ corresponds to the white noise component, aj is the rms intensity of the granules 
and Ti is the characteristic time scale of granulation. For the Gaussian term, the parameters 
-Pg, t'max, and (jg are the height, the central frequency, and the width of the power excess 
hump. 

Fig. [T] shows the power density spectrum of HD 186355, together with the fitted 
model using Eq. The three components of the background and the white noise were 
simultaneously fitted to a lightly smoothed power spectrum (Gaussian with FWHM of 
0.5 /iHz) outside the region where the power excess hump is seen. The value of z/max was 
obtained by fitting to a heavily smoothed power spectrum (Gaussian with FWHM of 3 
Ai/, where is estimated in Sect. 13. 2p . giving 106.5 ± 0.3 /iHz. Finally, the background 
and the white noise were subtracted from the power density spectrum, leaving only the 
oscillation signal (lower panel of Fig. [T]). 
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3.2. Individual Frequencies 

The background-corrected power spectrum in Fig. [1] shows the clear signature of 
solar-like oscillations: a regular series of peaks spaced by the large separation. We also see 
multiple peaks due to mixed / = 1 modes (see also Sect. IH Beck et al. 2011, Bedding et al. 
2011). The power spectrum is shown in echelle format in Fig. [21 both with and without 
smoothing. This diagram was made by dividing the power spectrum into six segments, each 
Az/ wide. We see that the peaks align vertically, allowing us to assign the / values indicated 
on Figs. [1] and [2j We do not see an obvious signature of rotational splitting, and the effect 
of stellar rotation is not considered in this paper. 

To ex tract the frequencies of individual oscillation modes, we used the software package 



Period04 (ILenz fc Breger 



20041 ). This uses iterative sinewave fitting, which does a good job 
of extracting frequencies in cases such as this, where the individual modes are unresolved 
or barely resolved. The red lines in Fig. |2] show the frequencies of 33 extracted peaks with 
SNR greater than 3). These are listed in Table [H together with their amplitudes and 
uncertainties. 

The uncertainties derived by Period04 are underestimates because they only consider 
the internal consistency of the parameters. We derived more realistic uncertainties (the 
second column in Table [1]) by means of Monte-Carlo simulations. The residual time-series 
{t, y) were obtained by subtracting the sum of multiple sine functions of frequencies 
(/i)i=i,2,...,n and the corresponding amplitudes Ai and phases 0i from the observed 
time-series {t, x) as 

y = x- ^ Aj sm{27r fit + (pi). (3) 

i=l,...,7i 

Then \y\ is regarded as the observational uncertainty in x. We constructed 100 simulated 
time-series z, which have the same residuals as the observed time-series. The 100 simulated 
time-series were fitted with the sum of multiple sine functions by taking (/j, Ai, 0j)j=i,2,...,n as 
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initial values according to least-squares algorithm. Hence 100 sets of new (/j, Ai, (/)j)i=i,2,...,n 
were obtained. The standard deviations of each parameter of (/«, ^i, 0i)«=i,2,...,n were then 
calculated, which were adopted as the uncertainty estimates of the parameters. 

We also calculated the mean large frequency separation Au by performing a linear fit 
to the five / = modes. Each data point was weighted by the uncertainty of the frequency 
listed in Table [1] Frequencies for / = modes are the most suitable ones for this calculation 
because they are not affected by the mixing with g modes. The slope of the fitted line gave 
the large separation to be Au = 9.37 ± 0.03 /iHz. 



4. Modelling 



The common way to estimate the fundamental properties is to compare calculated 
model parameters with the observational constraints. We employed the Yale rotating stellar 
evolution code (YREC; Demarque et al. 2008) for stellar evolution modelling computations, 
and the non-rad ial and non-adiabatic stellar pulsation programme JIG developed by 



Guentherl (119941 ) for frequency calculations. YREC can evolve our models up to the tip of 
the red giant branch, which is adequate for HD 186355. The input physics of the curr ent 
YRE C version (YREC7) includ ed the latest OPAL opaci ty tables fjlglesias fc Rogers 



996f). OPAL equat ion of state (IRogers fc Nayfonov 



( Angulo et al 



20021 ) and NACRE reaction rates 



19991 ). At low temperatures, opacities are obtained from 



Ferguson et al 



( 20051) . Convection is treated under the assumption of mixing length theory ( 



Bohm-Vitense 



19581 ). We did not take rotation, diffusion or convective overshoot into consideration in our 
calculation. 

There are several main inputs in YREC7 — mass, ami (to determine the mixing-length 
Imi = C(miHp, where Hp is pressure scale height), hydrogen abundance (X) and heavy- 
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element abundance [Z). The best models are searched among those grids after being 
compared with observational constraints. For our models, a^ni and X were fixed to the 
solar values of 1.8 and 0.72, respectively. The value of Z was varied within a certain range, 
usually from 0.005 to 0.025 with a step of 0.002, but it changes for models with different 
masses. 



Our initial estimate for the mass was made using scaling relations. iBrown et al.l (Il99l[ ) 
proposed a scaling relation that can be used to predict Umax by scaling from the solar case: 



M 

Me 



_R_ 

Rq 



'-eS 



eff,0 



-1/2 



(4) 



This relation gives a very good estimate for z/max for less evolved stars (iBedding fc Kjeldsen 



20031), while 



Stello et al 



(|2008[ ) ha ve shown that it holds also 



although with larger uncertainties, 
predict Au 

Av 



br stars on the giant branch. 



Kjeldsen fc Beddind (119951 ) give the scaling relation to 



M 



1/2 



R 



-3/2 



Az/0 \Mq) V^0, 
Knowing t'max, and Tefj, the stellar mass is estimated by: 



M 



Av 



-4 



cff 



3/2 



(5) 



(6) 



Using z/jnax and Az/ from Sec. [3] and Tgfj from KIC, the stellar mass is estimated as 1.41 ± 
0.14 Mq. Therefore, the initial masses of our models were chosen to be within the range of 
1.25 to 1.55 Mq with a step of 0.01 Mq. 

We looked for models for which the parameters are located inside the \-a error box 
confined by the uncertainties of observational results in the H-R diagram. For these sets of 
modelling parameters, we used a fine resolution for Z (in steps of 0.001) in order to find the 
best models. Some models with larger masses were also calculated, with a bigger mass step 
of 0.1 Mq, in an attempt to search for models in a large range, because the scaling relations 
and hence the estimated mass are not so reliable for the giant branch stars. 
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Fig. |3] shows several evolutionary tracks of models having different input parameters. 
The rectangle is the 1-a error box whose center corresponds to the observed stellar 
properties, from which we can see that HD 186355 is on the ascending giant branch, in the 
shell hydrogen-burning phase. Those models for which the parameters are within the error 
box and the mean large frequency separations are around 9.37 /iHz (within 0.03 /iHz) are 
indicated by dots. Different evolutionary tracks may pass through the same position in the 
H-R diagram by tuning the inputs. For example, a decrease of mass can be compensated by 
a decrease of hydrogen and heavy elements abundances to obtain the same position. Taking 
variations of the mixing-length into consideration, which move the tracks horizontally 
but have almost no influence on the luminosity, makes it even more complex to look for 
models. However, it is beyond the scope of this paper to consider the effects of varying the 
mixing-length and hydrogen abundance. We performed a minimization to find the best 
models. The definition of the function was based on two observed stellar parameters 
(luminosity and Tcs) and on the individual frequencies, as follows: 

2 /I r/r 1 r//r\2 ^ N / ,s 2 



2 _ i' Tcs-T^s \ f log L/Lq - log L'/Lq \ I ^ fuj-v: 
^ ' 150K J V 0.11 I N a. 



1=1 



where terms with primes are observed values, is the number of observed frequencies and 



(jj is the uncertainty of each frequency. We did not find it necessary to 



the model frequencies to correct for near-surface effects (IKjeldsen et al 



am 



20081) 



y an offset to 



We began by fitting only one mode of each degree in each order. For / = 1 we took the 
strongest peak in each order. The results based on these 17 observed frequencies are listed 
in Table [21 Although tracks for models with masses larger than 1.60 Mq also pass through 
the error box, their oscillation parameters differ greatly from the observations, which leads 
to bigger x^- From Table [2] we can see the best model has a mass of 1.43 Mq and Z of 
0.012. However, the difference between this and the 1.60 Mq model is very small, which 
makes the latter one a candidate for the upper limit of the stellar mass estimation. In order 
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to investigate this, we plot a series of echelle diagrams in Fig. H] to compare the theoretical 
and observed frequencies for each model shown in Table [2J Only those frequencies having 
minima of mode inertia (see Fig. E]) are shown, because those modes will have the highest 
amplitudes at the stellar surface. Good agreement is found for 1.43 M©, and models with 
larger masses do not reproduce the observed frequencies as well. However, the 1.60 Mq 
model is an exception which produces a rather good match to observed frequencies for / = 
and 2 modes. The location of the 1.60 Mq model is close to the center of error box in HR 
diagram which, combined with the relatively good fit to the / = and 2, leads to a smaller 
result than models with higher masses. After being compared with the theoretical modes 
in echelle diagrams, the degrees of observed modes (see Fig. [2]) are confirmed, including two 
modes with / = 3. 



As mentioned in Sec. 13. 2[ we found multiple oscillation peak s per order 



(see Table [B 



t hat w e interpret as mixed modes. As discussed by 



Bedding et al. 



Beck et al. 



or / = 1 



fl201lh and 



(I2OIII ). it is impossible to observe some mixed modes (g-dominated mixed 
modes) because they have very high inertias. However, other mixed modes act more 
like p modes (p-dominated mixed modes), having a lower inertia than th e g-dom ir iated 



mix ed modes and hence larger amplitude, which makes them observable. 



and 



Miglio et al. 



Tassouil mm 



(120081 ) have shown that pure g modes are equally spaced in period. 
P-dominated mixed modes are only approximately equally spaced in period. Measuring 
the period spacings for these observed mixed modes allows us to probe the cores of red 



giant stars. 



Beck et al. 



(120111 ) have detected mixed mo des in a red giant star with Kepler 



data and measured their period spacing. Subsequently, iBedding et al.l (1201 li ) have found 
a way to distinguish between hydrogen-burning and helium-burning red giants by using 
their different period spacings. They found that hydrogen-shell burning stars have observed 
period spacings mostly around 50 s, while stars with helium-burning cores have observed 
period spacings of about 100 to 300 s. The observed mean period spacing of HD 186355 is 
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58 ± 4 s which we m easured by means of power spectrum of the power spec trum method 



( Bedding et al 



201 ll ). This agrees with the value of 56 s for this star found by 



Bedding et al. 



( 120111 ). and confirms that HD 186355 is still in the shell hydrogen-burning phase. This 



value also agrees with our models. We note that measuring the period spacings may also 

provide a method to determine the size of the convective core of those helium-burning red 
I II ■ 



giants (iChristensen-Dalsgaardi i20 1 ll) . 



To make use of this extra information, we took frequencies of those / = 1 mixed modes 
with relatively low theoretical mode inertias into calculation. About 12 / = 1 modes were 
used for each model, giving a total of around 24 frequencies. These produced the values 
labelled Xi Table [2l Again the best model is the one with 1.43 Mq, but now the models 
with higher masses have large deviations between observed and theoretical oscillation 
frequencies. In particular, the 1.6 Mq model is ruled out after this calculation. We search 
for models with xl smaller than 30 and determine the mass to be 1.45 ± 0.05 Mq. 



5. Conclusion 

We have analysed the time series data sets of the star HD 186355 from Kepler to 
obtain its oscillation parameters. By using the scaling relations between Au, z/max and 
the stellar effective temperature Tgfj we estimated the stellar mass as 1.41 ± 0.14 Mq. In 
order to determine the stellar global properties more accurately, we computed a set of 
models to compare with the observational constraints. The best model was found having 
a mass of 1.43 Mq, which agrees with the scaling value, and Z of 0.012, and the stellar 
mass is constrained to be 1.45 ± 0.05 Mq. Furthermore, parameters such as age, effective 
temperature, luminosity and radius are also determined after comparison (see model with 
mass of 1.43 Mq in Table |2]). We also obtain the observed mean period spacing of / = 1 
modes with a value of 58 ± 4 s. From the modelled evolutionary track of HD 186355, we 
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know it is in the shell hydrog en-burning pha s e and on the ascending giant branch, which is 



consistent with the results of 
modes for red giants. 



Bedding et al. 



( I2OIII ) on the mean period spacings of mixed 
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Table 1. Frequencies extracted by Period04 



I Frequency Frcq. sigma Amplitude SNR 





(/.Hz) 


(mHz) 


(ppm) 




1 


81.5 


0.15 


0.37 


3.1 


2 


84.9 


0.31 


0.50 


4.2 





86.4 


0.56 


0.33 


3.0 


1 


90.7 


0.26 


0.79 


6.6 


1 


91.7 


0.25 


0.38 


3.1 


2 


94.1 


0.19 


0.61 


4.4 





95.4 


0.09 


0.61 


4.8 


1 


99.6 


0.11 


0.88 


6.8 


1 


100.2 


0.39 


0.97 


7.1 


1 


100.7 


0.30 


0.55 


4.2 


1 


101.3 


0.06 


0.67 


5.9 


1 


101.7 


0.22 


0.47 


4.5 


2 


103.6 


0.21 


0.75 


6.8 





104.7 


0.12 


0.96 


6.6 


3 


106.6 


0.66 


0.38 


3.2 


1 


107.7 


0.47 


0.54 


3.2 


1 


109.0 


0.86 


0.56 


3.7 


1 


109.3 


0.35 


0.67 


5.8 


1 


109.8 


0.58 


0.85 


6.7 


1 


110.2 


0.50 


0.55 


4.3 


1 


110.5 


0.29 


0.40 


3.0 


1 


110.8 


0.04 


0.38 


3.0 


2 


112.9 


0.21 


0.65 


5.7 





114.1 


0.25 


0.60 


4.6 
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Table 1 — Continued 



I 


Frequency 


Freq. sigma 


Amplitude 


SNR 




(/xHz) 


(/xHz) 


(ppm) 




3 


116.1 


0.25 


0.36 


3.0 


1 


118.8 


0.13 


0.45 


4.3 


1 


119.5 


0.46 


0.56 


4.6 


1 


120.4 


0.24 


0.36 


3.1 


2 


122.3 


0.12 


0.38 


3.2 





123.5 


0.66 


0.59 


4.5 


1 


127.0 


0.44 


0.33 


3.1 


1 


128.2 


0.11 


0.28 


3.0 


1 


128.8 


0.08 


0.33 


3.0 



Table 2. Modelling results: fundamental properties and for models with different masses. The expression of is 
given by Eq. ([7]), while Xi takes frequencies of mixed modes with low mode inertia into consideration. 



M/Mq 


Z 


Age 


Tcff 


L/Lq 


R/Rq 


log 5 


Au 




xl 






(Gyr) 


(K) 








(/xHz) 






1.35 


0.011 


3.26 


4886 


21.35 


6.45 


2.95 


9.37 


18.4 


44.5 


1.43 


0.012 


2.78 


4880 


21.98 


6.57 


2.96 


9.39 


2.1 


7.4 


1.45 


0.013 


2.76 


4859 


21.86 


6.61 


2.96 


9.37 


4.4 


29.8 


1.50 


0.011 


2.21 


4923 


23.56 


6.68 


2.96 


9.37 


2.8 


25.9 


1.60 


0.013 


1.85 


4895 


23.92 


6.81 


2.98 


9.37 


4.0 


60.5 


1.70 


0.019 


1.83 


4802 


23.06 


6.95 


2.98 


9.36 


26.1 


84.5 


1.80 


0.019 


1.49 


4825 


24.38 


7.07 


3.00 


9.37 


66.9 


97.0 


1.90 


0.018 


1.22 


4861 


26.16 


7.22 


3.00 


9.35 


36.6 


84.4 



Observational 
constraints 



4867 ± 150 24.0 ± 5.6 



9.37 ± 0.03 
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100.00 F 




Frequency (/i.Hz) 



Fig. 1. — Top panel: power density spectrum of the combined first five quarters of data (liglit- 
grey) and corresponding global model fit (black line). The dark-grey line is the smoothed 
(Gaussian with a FWHM of 0.5 /iHz) power density spectrum. The dotted line is the fitted 
background and the dashed line is the white noise. The red line shows the contribu tion from 



20081 ). Bottom 



stellar activity, the blue from granulation and the green from faculae (iKarofH I: 
panel: background corrected power density spectrum in the range of the stellar oscillations. 
It is clear that z/max is around 106 yuHz and that peaks are regularly spaced with a large 
spacing of 9.37 /iHz. Numbers are degrees of the modes. The multiple peaks correspond ing 



tn mnrlps nf / = 1 nrp> Kp>lipivp>rl +n Kp> mivfA mnrlpis rtRppV pi nl 901 1 • Rprlrlincr ft nl 901 1 ^ 
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Fig. 2. — The echelle diagram of the smoothed power density spectrum (dark, FWHM of 0.1 
/iHz) and the unsmoothed background corrected power density spectrum (grey, power divided 
by a factor of 2) divided into bins each Az/ wide. The red bars indicate 33 frequencies hsted 
Table [H The peaks for the same degree almost line up. The offset from perfect alignment 
is caused by the variation of the large frequency spacing with frequency. Multiple peaks for 
/ = 1 modes can be seen clearly. 
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Fig. 3. — Evolutionary tracks for eight models with different initial masses and heavy element 
abundances indicated by different colours listed in Table [2l The rectangle is the la error box 
for the observational constraints, and dots are models which have a large frequency spacing 
close to 9.37 /xHz. 
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Fig. 4. — Echelle diagrams for models plotted in Fig. |3l Squares, diamonds, triangles and 
circles are used for modes of degree Z = 0, 1, 2, and 3, respectively. Observed frequencies 
are indicated by filled triangles. For theoretical frequencies, only those with a corresponding 
observed mode are shown (except for / = modes) and symbol size indicates expected 
amplitude of each mode, which is scaled from mode inertia. For the observed ones, the 
size is scaled from the amplitude listed in Table [H Models with masses around 1.43 Mq 
reproduce the observed frequencies better than those with higher masses. The differences 
between individual frequencies are included in the calculation of xl- 
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Fig. 5. — Mode mass obtained by JIG9 for each degree (squares linked by solid line: / = 0, 
diamonds linked by dashed line: 1 = 1, triangles linked by dash dot line: / = 2, circles 
linked by dotted line: / = 3) versus frequency for the model with mass of 1.43 Mq that has 
the smallest x^- The dimensionless mode inertia is the ratio of mode mass to stellar mass. 
Those frequencies with minimal mode masses (hence minimal mode inertias) are plotted in 
the echelle diagram in Fig. |H 



